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ABSTRACT 

We present a result of age estimation for star clusters in M33. We obtain 
color-magnitude diagrams (CMDs) of resolved stars in 242 star clusters from 
the HST/WFPC2 images. We estimate ages of 100 star clusters among these, 
by fitting the Padova theoretical isochrones to the observational CMDs. Age 
distribution of the star clusters shows a dominant peak at log{t) ~ 7.8. Majority 
of star clusters are younger than log{t) = 9.0, while ten star clusters are older 
than log{t) ~ 9.0. There are few clusters younger than log{t) = 7 in this study, 
which is in contrast with the results based on the integrated photometry of star 
clusters in the previous studies. Radial distribution of the cluster ages shows 
that young to intermediate-age clusters are found from the center to the outer 
region, while old clusters are distributed farther from M33 center. We discuss 
briefly the implication of the results with regard to the formation of M33 cluster 
system. 

Subject headings: galaxies: star clusters — galaxies: individual (M33) — galaxies: 
photometry 



1. Introduction 

Ages of star clusters in a galaxy are important information for understanding the for- 
mation and evolution of star clusters and their host galaxies. In general we estimate ages of 
star clusters whose stars are resolved by fitting the color-magnitude diagrams (CMDs) of the 
stars with theoretical isochrones, and ages of unresolved clusters by fitting their integrated 
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photometry or spectra with population synthesis models. Star clusters in the Local Group 
galaxies are partially resolved into individual stars in the images obtained with the Hubble 
Space Telescope (HST) so that they play an important role in the study of extragalactic star 
clusters. 

There have been numerous studies that estimated the ages of star clusters in the Local 
Group galaxie s with the photometry of resolved stars in s tar clusters: for ex ample, LMC (e.g. 
Mackev et all J2006h and refere nces therein) and SMC Jciatt et al.lboosh. NGC 1 85. NGC 
205 Jsharina et allbooeh . M31 jRich et al.ll2"o05l ). and M33 JSaraiedini et al.l[200ol . boO?! ). 



M33 (NGC 598), Triangulum Galaxy, is a late-spiral galaxy (Scd), and is one of the 
three spiral galaxies in the L o cal Group. It is kii o wn to have numerous st ar clusters (see 



Sarajedini fc Mancone I (120071 ): iPark fc Lee I (120071): IZloczews. 
therein.). It is located at about 900 kpc from us (jCalleti et al. 



d et al.l (120081 ) and references 



20041 ) and has a large number 



of star clusters so that it offers us an excellent opportunity to study the ages of resolved star 
clusters. At the distance of M33, bright star clusters are easily resolved into stars and even 
the faint star clusters can be seen with a hint of resolved stars on the images taken with 
HST/Wide Field and Planetary Camera 2 (WFPC2). 



^There were several studies on the ages of star clusters in M33. IChandar et al.l (1999bl . 

200ll ) estimated the ages of M33 star clusters in the catalog of IChandar et al.l (Il999al . I2OO1I ) 
by comparing the integrated photometry of th e star clusters with a s imple stellar population 
(SSP) model in the color-color diagram. Later IChandar et al.l (|2002| ) supplemented their age 
estimates for more M33 star clusters with the spectroscopic data. They found that the age 
distribution of star clusters in the halo of M33 shows a large age spread, and little progression 
of abundance with age. They pointed out these properties are consistent with accretion origin 
for a large fraction of M33 halo but not with monolithic collapse scenario of the halo. From 
the comparison of globular cluster properties in nearby dwarf galaxies, they claimed M33 
halo cluster population might have formed by accretion of lower luminosity dwarf irregular 
galaxies and dwarf spheroidals rather the ones such as Magellanic cl ouds. The apparent lack 
of metallicity gradient among the ten M33 halo cluster samples of ISarajedini et al.l (120001 ) 
also further supported the accretion scenario. 



Ma et al.l estimated the a ges and me tallicity of the M33 star clusters in their series of 
papers jMa et all boOll . l2002al ffl. l2004al jbl) . They obtained the integrated spectral energy 



distribution (SED) of star clusters from the photometry using the Beijing-Arizona- Taipei- 
Connecticut (BATC) system with 13 narrow band filters. They fitted the SEDs with the SSP 
model to simultaneously determine the age and met allicity of star c l usters. They showed 
that their age estimations for star cluster samples from IChandar et al.l (Il999all200ll) are con- 



sistent with previous estimates. Using the age estimates of M33 clusters given by 



Ma et al 
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Sarajedini fc Mancone I (120071 ) showed that the number of clusters in the age distribution 
appears to dechne with age with no obvious breaks or abrupt changes, being fit well by a 
single power law. 



Sarajedini et al.l (120001 ) estimated ages of ten resolved globular clusters in M33 using the 
isochrone fitting technique to the deep CMDs of resolved stars derived from the HST/WFPC2 
images. They estimated the metallicity from the slope of red giant branch, and found that 
M33 star clusters show a significant second parameter effect. They concluded from this result 
that these star clusters might be se veral Gyrs younger th an the Galactic globular clusters 
if the second parameter is the age. IChandar et al.l (I20061) estimated the age of a globular 
cluster, C38 in the hst of [Christian &: Schommer I (Il982l ). to be ^2-5 Gyrs o ld, ut ilizing the 
integrated spectrum and the resolved member star CMD. I Sarajedini et al.l (120071 ) found 24 
star clusters from the deep images of two HST/Advanced Camera for Survey (ACS) fields 
in M33, and derived ages for 20 of them by fitting theoretical isochrones to the features on 
the CMDs of resolved stars including the main-sequence turnoffs (^MS TO). The derived ages 
of the star clusters range from log(t) = 7.5 to log(t) = 10.1. Recently IStonkute et al.l (120081 ) 
derived an age of t > 7 Gy r, for one extended clus t er at the outer edge of M33, using the 
CMD of the resolved stars. IZloczewski fc Kaluzny I (120091 ) also estimated ages of four M33 
star clusters by fitting isochrones to the member star CMD derived from HST/ACS images. 



Park &: Lee I (120071 ) compiled a comprehensive catalog of 242 star clusters in M33 identi- 



fied from the HST image data. Along with 32 new star clust ers detected from 24 HST/WFPC2 



fields, they combined several previo us M33 star catalogs (IBedin et al.ll2005l : IChandar et al. 



1999al . I2OOII : ISarajedini et al. 



20071 ) obtained with HST data only, thus to ensure that the 



resulting master catalog contains only the confirmed star clusters. This catalog offers a good 
sample of star clusters for age studies. 

As part of our project on studying star clusters in M33, we derived photometry of the 
resolved stars in and arou nd the star cl u sters from the HST/WFPC2 images that contain 
M33 star clusters hsted in iPark &: Lee I (120071 ). We estimated ages of 100 star clusters in 
M33 using the resulting CMDs of resolved stars. 

Section 2 describes the HST/WFPC2 data used in this study, and the data reduction 
procedures. Section 3 present the CMDs of M33 star clusters, and the estimates of cluster 
ages. In Section 4, we describe the age distribution of M33 star clusters thus derived, together 
with their integrated photometric parameters and their radial distributions from the center 
of M33. In Section 5 we discuss differences among the age estimations for M33 star clusters, 
and the implications of the results with regard to the formation of M33 cluster system. 
Finally, summary and conclusion are given in Section 6. 
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T hroughout th i s pap er, we will refer to the star clusters according to IDs in the cata- 
log of Park &: Lee I (2007 ). For cross-identifications, we will refer to the star clusters frona 
the c atalog of iBedin et al.l (120051 ) with B designation, and o nes from IChandar et al.l (Il999al . 
2001 ) with CBF. CS designation means t he star clusters fro m Christian fc Schommer J IQsil) , 



MKK SS designation means the ones from lMocheiska et al.l JlQOsh. MD the ones from 



Melnick fc D'Odorico 



(119781) ■ and SBGHS des i gnatio n means the star clusters from lSarajedini et al.l (120071 ). Note 
Sarajedini fc Mancone I (120071 ) compiled all pubhshed lists of M33 star clusters to build a 
comprehensive catalog which is being maintained on-lin^. For reference, we included their 
catalog IDs with SM designation. 



Data & Data Reduction 



We used HST/WFPC2 images of M33 fields available in the HST arch i ve. T hese 



HST/WFPC2 image set s include the 20 fields listed in Table 1 of lChandar et all (Il999al). 35 



f ields listed in Table 1 of lChandar et al.l (120011 ). and 24 fields listed in Table 1 in lPark fc Lee 
(120071 ). 



We derived the magnitudes of the point sources in the images using HSTphot (iDolphin 



2000al ) that is designed for the point spread function (PSF) fitting photometry of point 
sources in the WFPC2 images. HSTphot uses a library of TinyTim PSF for PSF fit- 
ting. Individual images were all coadded with utility program included in HSTphot. HST- 
phot produces ir istrum ental magnitudes as well as s tandard magnitud es calibrated with 
Holtzman et al.l ( 1l995l )'s relation and zeropoints from iDolphin I (l2C)00bl') after it au tomati- 
cally applies the correction for Charge Transfer Efficiency (CTE) ( iDolphin I l2000bl ) . Since 
HSTphot only calibrates the instrumental magnitudes of widely used filters, we used only 
the images obtained with FA36W{B), F555VI/(V), F606WiV), and F8UW{1). The resul- 
tant photometry data turn out to have small photometric error. Fig. [1] shows that resultant 
data have photometric errors smaller than 0.1 at y ~ 22 mag even for about 200 sec long 
image, and as for the data from exposures longer than 1000 sec, the errors get larger than 
0.1 only for V > 24. 



1 



http : //www . astro . uf 1 . edu/~{}ata/cgi-bin/m33_cluster_catalog/ index . cgi 
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3. Results 

3.1. Age Estimation 

Apparent sizes of the star clusters in the images are r ~ 1.0" in most cases, and as 
large as r ~ 3" for some cases with long exposures. We construct the CMDs of resolved 
stars in cluster-centric radii from 1" to 2.5" varying on the sizes of the clusters, which are 
listed in the fifth column of Table [H To minimize the field star contamination, the CMDs of 
field stars were compared with those of member stars around each star cluster. Field stars 
were selected in the annulus of which the inner radius is r ~ 3" and of the same area as the 
member star selection circle. 

Since HST/WFPC2 data used in this study were obtained from various observational 
projects, the depths of the images are not homogeneous. Thus, it happened that a star 
cluster observed on more than one field in our image sets resulted in different individual star 
detection on each appearance. In such a case, we used the CMD that has the largest number 
of detected stars for the analysis. 

Fig. [2] displays the CMDs of resolved stars in 100 star clusters thus derived. Some CMDs 
include only a few resolved stars. Nevertheless they show features that can be used for age 
estimation such as bright MS in young star clusters. The most reliable age determination of 
star cluster can be made with the location of the MSTO. We regarded the position of the 
brightest MS star as MSTO for young star clusters. However, even the deepest image set 
among these could not reveal the positions of the MSTO for old star clusters. In those cases, 
we had to rely on the RGB feature which is not sensitive to the age. Therefore, we estimated 
only the lower age limits of the star clusters by fitting the isochrones to the observed part 
of the RGB. 



We used the theoretical i sochrones in the Padova models (iGirardi et al.ll2002l ) to estimate 



the age of each star cluster. IChandar et al.l (jl999bl ) reported that metallicities of their 60 



star clusters span a large range from Z = 0.0002 to Z = 0.03, and that their young cluster 
samples can be well represented by the Z = 0.004 model, while there is a distinct turnoff 
to lower metallicity(Z = 0.001) around log(t) ~ 9.0. Accordingly, we used the model for 
Z = 0.004 for star clusters that turned out to have mainly blue or intermediate color stars 
and few red stars. As for the old star clusters whose detected members are mainly RGB 
stars, we used either model of Z = 0.004 ot Z = 0.001, whichever me t allicit y could fit the 



RGB better. We used the reddening values given in the iPark fc Lee I ( 120071 ) for isochrone 
fitting for each star cluster. 



We adopted as the distance to M33, 910 kpc ((m — M)q = 24.8) given by iKim et al. 
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( 120021 ) ■ Note there is a good summary of r e cent estima t es of the distance to M33 obtained in 
various ways given by (iGalleti et aLll2004j ). iKim et al.l ((2002) estimated the distan ce to M33 
using the /-band magnitude of the tip of red giant branch (TRGB) of field stars (iLee et al. 
19931 ). During this study, we visually checked that the /-band magnitude of TRGB for M33 
field stars in our images is indeed close to their estimation. Therefore we adopted their result 
as the distance to M33. 

Isochrones were fitted visually, and the upper and lower age limits of a star cluster were 
also estimated visually. Table [T] lists the ages of the M33 star clusters we derived, along with 
metallicities and the reddening values used. 



3.2. Age distribution 



Fig. El shows the age distribution for 100 clusters including the ten star clusters with 
lower age limit only, i n comparison with those in p revious studies: (a) Chandar sample for 
103 star clusters from lChandar et al.l (Il999bl . |2002| ) wh o used int e grated photo met ry for a ge 
estimation, (b) Ma sample for 226 star clusters from iMa et al.l J200ll . l2002al lblfl l2004al jbl) 
who also used in tegrated photometry fo r age estimation, and (c) Sarajedini sample for 20 
star clusters from lSarajedini et al.l (120071 ) who used the CMD of the resolved stars in clusters 
for age estimation. 

The age distribution of star clusters in this study shows a broad distribution with a 
peak at log{t) ~ 7.8, and a small number of star clusters older than log{t) ~ 9.0. The broad 
component can be fit well with a Gaussian with the peak age at log{t) = 7.83 ± 0.03 and 
the width of cr = 1.17 ± 0.11. Thus a majority of the clusters in this study are young to 
intermediate- age and a small number of the clusters are older than about one Gyr. It is 
noted that there is only one star cluster younger than log{t) = 7 in the our sample, while 
the Ma sample shows a large component of young clusters with a peak at log{t) = 7 and the 
Chandar sample shows a small component of young clusters with a peak at log{t) = 6.7. 

Due to the different depths of HST/WFPC2 images used in this study, it is not easy 
to estimate the completeness of the age sample. However, we note that the isochrone fitting 
method used in this study has a bias such that bright, large sparse star clusters are likely to 
be included more than faint compact ones in the final age sample. 
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3.3. Integrated color-magnitude diagram 
Fig, m d i splays the CMD of integrated photometry for 242 star clusters in M33 from 



Park fc Lee I (120071 ). along with the information of age derived in this study. In general, the 
age and integrated magnitude/color show consistent behaviors on the CMD. Old star clusters 
with log{t) > 9.0 have colors —0.5 < {B — V)o < 1.1 and magnitudes —6 < My < —9, which 
are similar to those of the Milky Way globular clusters (MWGCs), 0.6 < {B — V)o < 1.0 and 
— 10.0 < Aly < —6.0. Therefore, these old star clusters can be regarded as being analogous 
to the MWGCs. 

One notable star cluster is cluster ID 85 (= MKKSS-27, SM-261), with My = -9.114 
and [B — V)o = 0.578. It is the brig htest in the sani p le. T hese values are typical magnitude 



and color of an old globular cluster. IChandar et al.l (120021 ) estimated its age to be log{t) 



8.6 ± 0.2 from integrated photometry, or 2-5 Gyrs from spectroscopy. We obtained PSF 
photometry of ten stars around this star cluster. Summing the magnitudes of these ten stars, 
we derive the total magnitude and color V = 16.051 ± 0.002 and {B - V) = 0.657 ± 0.004. 
These values are very simila r to those from ap erture photometry, V = 16.017 ± 0.006 and 



B -V) = 0.678 ± 0.004 (IPark fc Lee 1 120071 ). Therefore we regarded PSF photometry 



recovered most of bright member stars in this cluster. The HST/WFPC2 images of the field 
including this cluster and the CMD of this cluster show that there are two very bright red 
stars. The spatial locations of these red stars in the cluster indicate that they are probably 
member stars, and their locations on the CMD show that they are red supergiant stars. The 
fluxes of these stars contribute to most of the integrated light of the star cluster. If we remove 
the light contribution of these two red supergiant stars, we get the integrated photometry 
of y = 18.625 ± 0.015, and (B-V) = -0.004 ± 0.019. This shows that this cluster is much 
younger than the age derived from the case including the two stars. We estimated the age of 
this star cluster to be log(t) = 7.0, which is much younger than previous estimations. This 
case clearly shows the power of age estimation with resolved star CMD over those based on 
integrated light. 

Intermediate-age star clusters (7.8 < log{t) < 9.0) occupy the intermediate color region 
in the CMD. It is interesting to note that their integrated magnitude range is quite fainter 
than those of younger star clusters or of older red star clusters. On the other hand, star 
clusters younger than log{t) ~ 8.0 have bright young stars whose fluxes contribute the 
majority of total star cluster light. These stars are bright main sequence stars or supergiants. 
Consequently, young star clusters are located in the bright, blue regions on integrated CMD. 

Since the ages in this study were derived in an independent way from using the integrated 
color of the star clusters, more insight to star cluster properties can be obtained by comparing 
the ages and the integrated colors of the star clusters. Fig. displays the integrated color 
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versus age diagram. We plotted also the theo retical relations derived from single stellar 
population models of ISruzual fc Chariot I J2003h . for Z = 0.02, 0.004 and 0.0004. The figure 
shows that the integrated colors of the M33 clusters increase slowly as the age increases, 
approximately consistent with the theoretical curves. {V — I) and {B — /) colors of the 
young clusters show a large scatter for young star clusters, while their {B — V) colors show 
a small scatter around the theoretical curves. It is suspected that this may be due to the 
large uncertainty in the J-band photometry of the clusters. It is needed to obtain better 
integrated photometry in the J-band. 



3.4. Radial distribution 



With 100 star cluster samples of known ages in hand, we investigated their radial dis- 
tribution around M33. With consideration of their integrated {B — V)o colors and the age 
distribution shown in Fig. [31 we divided the star clusters into three age groups: (1) young 
star clusters with log{t) < 7.8, (2) intermediate-age star clusters with 7.8 < log(t) < 9.0, 
and (3) old star clusters with log{t) > 9.0. The boundary value of log{t) = 7.8 divides the 
young population in half seen in Fig. [3l and this age corresponds to the color {B — V)o ~ 0.2. 

Fig. shows the number of the clusters as a function of the p rojected radial distan ce 
from the center of M33 (a(2000) =1:33:50.9, 5(2000) =30:39:35.8, Jskrutskie et allbooeh l 
Although it is not easy to investigate with small number of samples, there exists a tendency 
that older star clusters are distributed in a wider area than the young star clusters in M33. 
Only young star clusters are found in the region close to the M33 center, whereas only 
intermediate-age a nd old star clusters are f ound in the outside region. This is consistent 
with the finding by ISarajedini fc Mancone I (120071 ) who used the age estimates derived from 
integrated photometry. 



4. Discussions 



4.1. Comparison with other studies 



Sarajedini et al.l (120071 ) estimated the ages of 20 star clusters by fitting isochrones to the 



stellar photometry data derived from HST/ACS images. ISarajedini et al.l (120071 ) 's photom- 
etry is deep enough to reach the MSTO for most of their target s tar clusters. The r e are six 
star clusters with age estimates common between this study and ISarajedini et al.l (120071 ) 's: 
Clusters ID 120 (CBF-110, SM-103), 130 (CBF-108, SM-127), 131(CBF-106, SM-130), 228 
(CBF-70, SM-402), 231 (CBF-79, SM-409), and 232 (CBF-78, SM-410). Fig. [7] displays the 
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compariso n of two age estimates fo r these star clusters. The age estimates show a good agree- 
ment with ISarajedini et aLl (120071 ) 's for four clusters. We derived only the lower age limits 



for tw o clusters: log{t) > 9.4 for ID 228 and log(t) > 9.2 for ID 131, while ISarajedini et al. 
J2007h derived log{t) = 10.1 ± 0.05 and 9.8 ± 0.05, respectively. 



Age estim a tion of M33 star clusters with integrated photometry was carried out by 



Chandar et al.l (jl999bl . |2002| ) for most of their 154 star cluster samples. IChandar et al. 



(2002) obtained ages for additiona l 14 star cluste r s frorn the catalog of iMochejska et al. 



(119981 ) with the same method as in I Chandar et al.l (Il999bl ). in addition to revising the age 
estimates for some clusters with spectroscopy. There are 58 star clusters with age estimates 
common between theirs and this study. Fig. El^b) shows that both age estimations agree 
within the error of measurements for many of star clusters. But our age estimates for star 
clusters of log{t) = 8.0 ~ 9.0 are younger than theirs in general. On the other hand, there are 
considerable differences for young clusters such that our estimates are mostly log(t) > 7.0, 
while theirs are log{t) < 7.0. 

We carefully ex a mined the st ar clusters for which our age estimates are much different 
from I Chandar et al.l (jl999bl . |2002| ). Fig. [H] shows the CMDs of ten young star clusters that 
have large age differences. To compare the age estimations, we plotted the Padova isochrones 
for both age estimates on each CMD. y-band grayscale images of the star clusters are also 
displayed to show the spatial distribution of detected stars. It is seen that there are some 
very bright stars that were not detected by HSTphot in clusters ID 143 and 232 (SM-168 
and SM-410, respectively). They might have led to younger age estimation of star clusters 
had they been detected and measured. There are some undetected stars in cluster ID 170 
(SM-220), but they look like faint stars. Except for above three star clusters, it seems that 
HSTphot detected all bright stars in the region of each star cluster. Their resolved star CMDs 
are more consistent with isochrones for our age estimates than those for their estimates. In 
short, seven out of ten star clusters show features that are more consistent with our age 
estimations th an with theirs. Theref ore, we argue that our age estimations are more reliable 
than those by I Chandar et al.l ((2002) for most of star clusters with a large discrepancies in 
the age estimations. 

There is also considerable amou nt of age discrepancie s for some older star clusters: our 
age estimates are log(t) < 9.0, while IChandar et al.l (120021 ) estimated their ages to be much 
older. The CMDs and V grayscale images of these star clusters are shown in Fig. [91 Clusters 
ID 49, 53, 111, and 231 (SM-45, SM-95, SM-110, and SM-409, respectively). turn out to have 
only red star sequences in their resolved star CMDs. It is difficult to derive reliable ages for 
these clusters just with RGB features. 



But for other star clusters, we noted that there are some blue stars whose colors {{V — I) 
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or {B — V)) are ~ 0.0, with their magnitude errors smaller than 0.2. Most of these stars are 
located in the very central region of the star clusters. Therefore, we regard the se as blue main 
sequen ce member stars. If they are indeed the member stars of the cluster, 
2002!) might have overestimated the ages of these star clusters. Note that 



Chandar et al. 



Chandar et al. 



( 1l999bl ) estimated from the integrated photometry the a.ges{log{t)) of clusters ID 35, 59, 169 



and 219 (SM-68, SM-146, SM-219, and SM-377, respectively) to be 8.3 ± 0.30, 7.7 ± 0.20, 
8.4 ± 0.20, and 7.8 ± 0.30, respectively, which are very close to our estimates: 8.3 ± 0.2, 
7.8 ± 0.3, 8.0 ± 0.2, and 8.1 ± 0.2. 

We also compare d our age estimates of the star clu sters with the results obtained in 
a series of papers by iMa et al.l J200ll . l2002al ffl l2004al Jbl). as shown in Fig. [7](c). There 
are 71 star clusters common between both studies. For star clusters older than log{t) ~ 
8.0, our age estimates are in general younger than their estimates except for a few star 
clusters. As for the star clusters of log{t) < 8.0, it turned out that our ag e estimates are 



much o l der th an their estimates, basically the same tendency as those with IChandar et al. 



(Il999bl. l2002r) with larger d ifferences. We note they obtained (i? — V) colors to be bluer 
than I Chandar et al.l (1l999a[ l for the same clusters (iMa et al.l 1200 ll ) . It seems that large 
seeing condition (~ 2.0") and large aperture (r = 6.8") used in their photometry caused 
light contamination from neighboring b lue disk stars i n M33, thus leading to younger age 
estimation even compared wit h those of IChandar et al.l . This can be seen in Fig. Wid) where 
age estimates by lMa et al.l and lChandar et al.l are compared. Although many of age estimates 
are in good agreement between them over a wide age range, quite a number of star clusters 



show large differen ces. For al most all of such star clusters, IChandar et al 



s estimates are 



older than those by iMa et al. 



Sarajedini fc Mancone I (120071 ) compiled the data in the eight previous catalogs of M33 



star cluster and built a new catalog of 451 star cluster candidates in M33, of which 255 are 
sta r clusters confirmed with HST images. Since their age estimation is based on the results 
by iMa et al.l . we do not show the result here but refer to Fig. Wic) for the age comparison 
between ours and theirs. 

These comparisons indicate that the estimated ages for some star clusters in this study 
have uncertainties due to the possible inclusion of non-member M33 field stars. This is true 
to both integrated photometry or PSF photometry of resolved stars. However, the effect of 
non member star selection can be reduced in resolved star photometry since the contribution 
of individual stars can be evaluated. For example, CMD of stars near star cluster center, 
and those of control field are compared to minimize the field star contamination. Although 
PSF photometry could not detect some bright stars in the cluster region, we argue that our 
age estimations for star clusters based on the isochrone fitting are more accurate than other 
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age estimations based on the integrated light of the star clusters. 



The resultant age distribution shows that there is only one star cluster which is younger 
than log{t) ~ 7.0, as shown in Fig. [3](a). It might be that there are indeed few young star 
clusters in M33. However, M33 is a late spiral galaxy, in which signatures for active star 
formation are clearly seen in a number of giant HII regions such as NGC 604 or NGC 595 
on its disk. Far infrared observation indicates M33 is more actively f orming stars (with sta r 
formation rate ~ 0.3 — 0.7MQ/year) than either our Galaxy or M31 (IHippelein et al.ll2003l ). 
Therefore it is unlikely that M33 lacks very young star clusters. I n fact, the presence of very 
young clusters that are 3-5 Myrs old was reported in NGC 604 (IHunter et al.lll996l ). Such 
young clusters normally are embedded in nebulous environments, which makes it difficult to 
be detected with visual investigation. 



4.2. Formation of old star clusters 



The interaction between galaxies is known to often trigger the formation of stars and 
star clusters. For example, a large population of young bright star clusters in the interacting 
galaxy pair of NGC 4038/4039 (Antennae galaxies) a re believed to be inv olved with the 
interaction of their host galaxies (IWhitmore et al.lll999l ) . iLee &: Lee I ( 120051 ) suggested that 
the intermediate-color clust ers in the interacting galaxy NGC 1487 were probably formed 
during the merging process. iLee et al.l (120051 ) found out that the age distribution of the star 
clusters in NGC 5194 shows a broad peak that is consistent with the crossing times of the 
companion galaxy NGC 5195 through NGC 5194 disk. These results indicate that the age 
distribution is consistent with the scenario that the interaction between two galaxies directly 
caused an increase in cluster formation rate. 

In the case of M 33, the spatial distribution of halo stars shows no sign of disturbance 
( iFerguson et al.ll2006l). However, H I surveys around M33 revealed an HI bridge between M31 
and M33 ( Braun &: Thilker 20041) . and the outer HI g as disk of M33 is found to be skewed 



towards the M31 direction (jCorbelli &: Schneider 1119971 ). From the model simulation. iBekki 



( 120081 ) suggested the tidal interaction between M31 and M33 about 4-8 Gyrs ago can explain 
the HI bridge and the observed HI warp in M33. There are ten old star clusters with lower 
age limits of log{t) ~ 9.0 in our sample. While we cannot estimate their exact ages with 



out data, previous studies on some of these clusters suggested tha t they are sever a. 
younger than Galactic g 



obula r clusters (see ISarajedini et al.l (120001 ): IChandar et al.l (120061 ): 



Gyr s 



Zloczewski fc Kaluzny I (120091 )). If the peak for the age distribution of these old star clusters 



are found to be similar to the epoch of the interaction, it might be possible that the galaxy 
interaction with M31 had triggered the formation of old star clusters in M33. We note that 



Huxor et al.l (120091 ) suggested the possibility that the tidal interaction between M31 and 
M33 have affected the spatial distribution of M33 globular clusters at large galactocentric 
distances. 

However, the small number of old star cluster samples prevents any definite conclusion 
on the picture of their formation. Further studies with more sample with good age estimation 
are needed to investigate in detail the formation of M33 star cluster system. 



5. Summary and Conclusion 



We present BVI photometry of the resolved stars in star clusters on the HST/WFPC2 
images of 79 fields in M33. Most star clusters in M33 are resolved into individual stars on 
the HST/WFPC2 images. We estimate the ages of 100 star clusters matching the CMDs of 
resolved stars with Padova theoretical isochrones. 

The age estimates for the star clusters in this study show a broad distribution from 
log{t) ~ 7.0 to log{t) ~ 9.0 with a peak at log{t) = 7.82 ± 0.03. There are also a small 
number old star clusters with lower age limits of log{t) ~ 9.0. It is found out that there are 
few clusters younger than log{t) = 7.0 in this study. 

Previous estimates Jchandar et~aDll999bl . I2OO2I : iMa et al.l[200ll . l2002al jblfl l2004al lb[) for 

some young star clusters show a difference from our estimates. Close inspection of the 
CMDs and images of each star cluster reveals that previous estimations with integrated light 
underestimated the ages for about 80% of young star clusters that have significant differences 
among the age estimations. 

Age distribution along the galactocentric distance reveals that there is a hint that blue 
young star clusters show a stronger concentration toward the M33 center than older star 
clusters. The interaction between M31 and M33 might affect the formation and spatial 
distribution of old star clusters in M33. 
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Fig. 1. — Photometric errors versus V magnitudes for a representative set of star clusters in 
this study. Upper four panels show the photometry with exposures longer than 1000 seconds 
while lower four panels show the results with exposures shorter than 1000 seconds. 
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Fig. 2. — V — {B — V) CMDs of the resolved stars in M33 star clusters. Filled circles 
denote the stars found within the member selection radius of each star cluster, and crosses 
denote the field stars around the star clusters. Solid lines represent the Padova theoretical 
isochrones fit to the data. ID and derived age of each star cluster are shown at the lower 
right corner and upper left corner of each plot, respectively. 
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Fig. 2.— V ~ (B-V) CMDs continued. 
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Fig. 2. — V — {V — I) CMDs of the resolved stars in M33 star clusters. 
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Fig. 7. — Comparison of a g e estimates in this stud y with tho s e from fa) ISaraiedini et al. 



(120071), (b) IChandar et al.l fllQQQbl . 120021), and (c) 



Ma et al 



Panel (d) shows t he comparison between IChandar et al.l (Il999bl . 



2002a 
limits. 



a 



2004 



3E) 



(12001 



2002a 



20021 ) and 



iraiedmi et al. 
blU l2004aljbh . 



Ma et al. 



(12001 



Open circles with arrows represent the star clusters with lower age 
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Fig. 8. — CMDs and y-ba nd grayscale im a ges of the star clusters for which our age estimates 
are larger than those by IChandar et al.l (120021 ). Intensity scales for grayscale maps are 
intentionally stretched so that the bright sources in each cluster can be clearly distinguished 
at the expense of faint sources. A large circle in each grayscale image denotes the member 
selection radius listed in Table [TJ Solid lii ies and dashed lines in the CMDs represent the 
Padova isochrones for our age estimates and lChandar et al.l (120021 ) 's. respectively. The points 
in each CMD represent the stars marked with open circles on the corresponding grayscale 
image. 



-27- 




Fi g. 9. — Same as Fig, p but for star clusters that were estimated to be younger by us than 
by lChandar et aL (|2002[ l. 
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Table 1. Age estimates for M33 star clusters 



ID^ 


log(t) 


Z 


E{B - y)^ 


r(arcsec) 


Cross identification'' 


101 




0.004 


0.15 ±0.05 


1.0 


CS-U70, SM-342 


102 




0.004 


0.20 ±0.05 


1.2 


CS-C33, MD-36, MKKSS-45, SM-351 


104 


o 0+0-2 


0.004 


0.10 ±0.05 


1.0 


CBF-154, CS-H19, SM-355 


108 


' -0-0.4 


0.004 


0.15 ±0.05 


1.0 


CBF-40, SM-35 


109 


o c+0.3 

o-o-o.e 


0.004 


0.15 ±0.05 


1.5 


CBF-86, CS-U140, SM-40 


111 


D 7+0.3 
^•'-0.2 


0.004 


0.20 ±0.05 


1.0 


CBF-53, SM-45 


112 


7 O+0.3 
' •'J-0.4 


0.004 


0.15 ±0.05 


1.0 


CBF-88, SM-47 


113 


> 9.4 


0.004 


0.10 ±0.05 


1.2 


CBF-54, CS-U137, MD-8, SM-49 


118 


8 

°-0-0.1 


0.004 


0.10 ±0.05 


1.0 


CBF-111, SM-100 


119 


> 9.1 


0.004 


0.10 ±0.05 


1.5 


CBF-114, CS-C38, MD-11, SM-102 



''Taken from lPark fc Lee I (|2007h 



^B identifications are from Bedin et al 



( 2005 ); CBF identifications from IChandar et al 



()l999al 12001); CS identificat i on fro m IChristian fc Schommer I (119821): MKKSS identifi- 
cation s from IMocheiska et all (119981); MP identificat i ons a re from iMelnick fc D'Odorico 
( 1978); SBGHS identifi c ation s from ISaraiedini et al.l (|2007l ): and SM identifications from 
Saraiedini fc Manconel |2007t) 



Note. — The complete version of this table is in the electronic edition of the Journal. 
The printed edition contains only a sample. 



